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Abstract 

We critically examine the definition of dijet cross sections at large rapidity intervals in 
hadron-hadron collisions, taking proper account of the various cuts applied in a realistic ex- 
■ perimental setup. We argue that the dependence of the cross section on the precise definition 

of the parton momentum fractions and the presence of an upper bound on the momentum 
transfer cannot be neglected, and we provide the relevant modifications to the analytical 
formulae by Mueller and Navelet. We also point out that the choice of equal transverse 
momentum cuts on the tagging jets can spoil the possibility of a clean extraction of signals 
of BFKL physics. 



1 Introduction 



Long ago Mueller and Navelet suggested Q to look for evidence of Balitsky-Fadin-Kuraev- 
Lipatov (BFKL) evolution g |, Q by measuring dijet cross section at hadron colliders as a 
function of the hadronic centre-of-mass energy y/S, at fixed momentum fractions x a j, of the 
incoming partons. This is equivalent to measuring the rates as a function of the rapidity interval 
Ay = \y a — Ub\ between the jets. In fact, at large enough rapidities, the rapidity interval is well 
approximated by the expression Ay ~ ln(s/|i|), where s = x a XbS and \t\ ~ k a ±kb±, with k a< b± 
being the moduli of the transverse momenta (i.e., the transverse energies) of the two jets. Thus, 
since the cross section tends to peak at the smallest available transverse energies, Ay grows as 
In S at fixed x a ^. 

It is clear that the measurement proposed by Mueller and Navelet is not feasible at a collider 
run at a fixed energy; on the other hand, to look for the BFKL-driven rise of the parton cross 
section directly in the dijet production rate da /Ay as a function of Ay is difficult due to the 
steep fall-off of the parton densities || |J. This led some of us to propose the study of less 
inclusive observables. In particular, the angular correlation between the tagged jets, which at 
leading order (LO) are back-to-back, is smeared by gluon radiation and by hadronization. Part 
of the gluon radiation originates from the mechanism responsible for BFKL effects, namely from 
gluon radiation in the rapidity interval between the jets. Accordingly, the transverse momentum 
imbalance || [7| and the azimuthal angle decorrelation || |(| ||, [| pi]] have been proposed as 
observables sensitive to BFKL effects. The azimuthal angle decorrelation has indeed been studied 
by the DO Collaboration at the Tevatron Collider JO] . As expected, a NLO partonic Monte Carlo 
generator, JETRAD 12, predicts too little decorrelation. However, the BFKL formalism 



predicts a much stronger decorrelation than that observed in the data. In fact, the data are well 



described by the HERWIG Monte Carlo generator []14|, [L5|, |iq| , which dresses the basic 2 — > 2 



parton scattering with parton showers and hadronization. This hints at a description of the 



azimuthal angle decorrelation in terms of a standard Sudakov resummation [17]. It therefore 
appears that, in the presence of Sudakov logarithms, it is quite difficult to cleanly extract the 
presence of BFKL logarithms from this observable, not least because the latter are expected to 
be smaller than the former in the energy range explored at present. It thus comes as no surprise 
that DO Collaboration [11] find no strong evidence of BFKL effects in their data. 



Recently, the DO Collaboration [18| has revisited the original Mueller-Navelet proposal, and 
has measured the ratio 

R = a{ ^} (1.1) 

of dijet cross sections obtained at two different centre-of-mass energies, \J~S^ = 1800 GeV and 
\JS B = 630 GeV. The dijet events have been selected by tagging the most forward/backward 
jets in the event, and the cross section is measured as a function of the momentum transfer, 
defined as Q 2 = k a ±kb±, and of the quantities 

X! = cosh(Ay/2) , x 2 = ^e~S cosh(Ay/2) , (1.2) 

with y = (y a + yfe)/2, Ay = y a — y^ > 0, and y a (yb) are the rapidities of the most forward 
(backward) jet. The dimensionless quantities x\ and x 2 are reconstructed from the tagged jets 
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x\ range 


%2 range 


bin 1 


0.06-0.10 


0.18-0.22 


bin 2 


0.10-0.14 


0.14-0.18 


bin 3 


0.10-0.14 


0.18-0.22 


bin 4 


0.14-0.18 


0.14-0.18 


bin 5 


0.14-0.18 


0.18-0.22 


bin 6 


0.18-0.22 


0.18-0.22 



Table 1: (x\-,%i) bins, with the upper bound of range in x± not larger than the upper 
bound of the range in x-i . 



using Eq. ( |1.2[ ), irrespective of the number of additional jets in the final state. In leading- 
order kinematics, for which only two (back-to-back) jets are present in the final state, we have 
x\ = x a and X2 = Xb, the momentum fractions of the incoming partons. Higher-order corrections 
entail that these equalities no longer hold; however, x\ = x a and X2 = Xb are still reasonable 
approximations (unless one goes too close to the borders of phase space). This implies that 
when the ratio in Eq. ( |1 . 1|) is computed at fixed x\ and X2, the contributions due to the parton 
densities cancel to a large extent, thus giving the possibility of studying BFKL effects without 
any contamination from long-distance phenomena. 



In the analysis performed by DO [18|, jets have been selected by requiring k a ,b± > 20 GeV, 
\lla,b\ < 3, and Ay > 2. The cross section was measured for ten (2:1,2:2) bins, of which we list 
m Table [l] the six with the upper bound of the range in x\ smaller than or equal to the upper 
bound of the range in X2 (the others may be obtained by interchanging x\ <-> X2). Finally, a cut 
on the momentum transfer, 400 < Q 2 < 1000 GeV 2 , has been imposed. With larger statistics, 
a binning in Q 2 would also be possible. These cuts select dijet events at large rapidity intervals. 
To have a crude estimate of the typical Ay values involved, we observe that, in a given (x\,X2) 
bin, the data accumulate at the minimum x\ and X2 in order to maximise the parton luminosity, 
and at minimum k± in order to maximise the partonic cross section. We can then use the 
LO kinematics (|1.2|) to obtain the effective rapidity interval. For instance, in bin 5 we find 
Ay A ~ 5.3 at -^5^=1800 GeV, and Ay B ~ 3.1 at a/S^~ = 630 GeV. In addition, we see that in 
the large-Ay limit, Ay A — ► Ay B + ln(S A /S B ). 

The data collected by DO are compared to BFKL predictions as given by Mueller and Navelet, 
and an effective 'BFKL intercept' is then extracted (see Eq. fl2.16| ) below). However, we argue 
in this paper that the different reconstruction of the x's used by DO as compared to the original 
Mueller-Navelet analysis (see Eqs. ( |1.2| ) and ( 2.11| )), and some of the acceptance cuts imposed in 
the experimental analysis, like the introduction of an upper bound on the momentum transfer 
Q 2 , actually spoil the correctness of this procedure, and require modifications of the Mueller- 
Navelet formulae. These modifications are subleading from the standpoint of the BFKL theory, 
however they have an impact on the extraction of the BFKL intercept at subasymptotic energies. 
Furthermore, the fact that dijet events are selected by means of transverse momentum cuts which 
are the same for the two tagged jets poses additional problems: large logarithms of (non-BFKL) 



2 



perturbative origin enter the cross section, and thus the ratio of Eq. (1.1) is affected by the same 
kind of problems as the azimuthal decorrelation. 

In this study we address the quantitative importance of these issues on the DO analysis using 
a combination of analytic and numerical techniques and several different theoretical approxi- 
mations for dijet production. The paper is organized as follows: in Section ^, we study the 
impact of a different definition of the x's and of the upper bound on Q 2 in the framework of 
the naive BFKL equation, and compare the result with the standard Mueller-Navelet analysis. 
In Section ||, our BFKL predictions are refined by considering running-coupling effects, and 
implementing energy-momentum conservation; this is done by using Monte Carlo techniques. In 
Section |4] the problem of equal transverse momentum cuts is investigated; in Subsection 4.1 



wc 



use a fixed-order perturbative computation (accurate to NLO), and in Subsection 4.2 we return 
to the BFKL formalism. Finally, Section || reports on our conclusions. 



2 Dijet production at fixed x's in the naive BFKL approach 

In large-rapidity dijet production, the rapidity interval is related to the kinematic invariants 
through Ay ~ ln(s/|i|), with s the squared parton centre-of-mass energy and \t\ of the order of 
the squared jet transverse energy. Thus, in computing the dijet production rate we encounter 
large logarithms, which arise in a perturbative calculation at each order in the coupling constant 
ot, s . The large logarithms can be resummed in leading-logarithmic (LL) approximation through 
the BFKL equation. 

In the high-energy limit, s 3> the BFKL theory assumes that any scattering process is 
dominated by gluon exchange in the crossed channel. This constitutes the leading (LL) term of 
the BFKL resummation. The corresponding QCD amplitude factorizes into a gauge-invariant 
effective amplitude formed by two scattering centres, the LO impact factors, connected by 
the gluon exchanged in the crossed channel. The BFKL equation then resums the universal 
LL corrections, of C(a™ ln n (s/|t|)), to the gluon exchange in the crossed channel. These are 
obtained in the limit of strong rapidity ordering of the emitted gluon radiation, i.e. for n + 2 
gluons produced in the scattering, 

Va >?/l >2/2 » ••• »2/n-l »l/n » 2/6 • (2-1) 

In dijet production at large rapidity intervals, the crossed-channel gluon dominance [] ensures 
that the functional form of the QCD amplitudes for gluon-gluon, gluon-quark or quark-quark 
scattering at LO is the same; they differ only by the colour strength in the parton-production 
vertices. We can then write the cross section in the following factorized form |l]] 

= / dkl ± dklJ^ X l,l)fM4^l)^^, (2-2) 

*The crossed-channel gluon dominance is also used as a diagnostic tool for discriminating between different 
dynamical models for parton scattering. In the measurement of dijet angular distributions, models which feature 
gluon exchange in the crossed channel, like QCD, predict a characteristic sin -4 ((9*/2) dijet angular distribution 
|ic| , |2l| l , while models featuring contact-term interactions, which do not have gluon exchange in the crossed channel, 
predict a flattening of the dijet angular distribution at large s/\i\ Eg| . 
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where fx F is the factorization scale, and x®, x® are the parton momentum fractions (1.2) in the 
high-energy limit, 

x° a = *^ e *« x° b = ^e-w , (2.3) 

and the effective parton densities are |24]] 

(2.4) 



f cS (x,fi 2 F ) = G(x,nl) + g [Qf&'Vl) +Qf(x,t4)] , 



f 



where the sum runs over the quark flavours. In the high-energy limit, the azimuthally-averaged 
gluon-gluon scattering cross section becomes M 



da, 



.</.</ 



dkl ± dk 2 ± 



k 2 



f(q a ±,qb±,Ay) 



k 2 
K b± 



(2.5) 



where Ca = N c = 3, and q a ±, q b ± are the momenta transferred in the t-channel, with q a ± = —k a ±_ 
and %± = kb±. The quantities in square brackets are the LO impact factors for jet production. 
The function / (q a ±, </&_l, Ay) is the azimuthally-averaged Green's function associated with the 
gluon exchanged in the crossed channel, and is obtained by integrating out the azimuthal angles 
in the solution of the BFKL equation, 



with 



f(q a ±,qb±,&y) 



1 



dv e' 



2 \ w 
j{v)Ay ( q a ± 

qL 



-2a s [Re^(l/2 + ^) - V(l)] 
= A-Bu 2 + 0{v A ), 

with ip the digamma function, and 

a s = a s N c /n, A = 4a s ln2, B = 14C(3)a s 



(2.6) 



(2.7) 
(2.8) 



(2.9) 



2.1 The standard Mueller-Navelet analysis 



In order to elucidate how the DO Collaboration |18| evaluates the effective BFKL intercept, we 
follow the original Mueller-Navelet approach p]], i.e. we substitute Eqs. ( |2.5| ) and ( |2.6| ) into 
Eq. ( |2.2[ ) and integrate it over the transverse momenta fc a j_ and k b ± above a threshold E±, at 
fixed coupling a s and fixed x®, x®. The rapidity interval Ay = \y a — yb\ i n Eq. ( |2.6| ) is determined 
from the x's (|2.3|), 

Ay = l Q pi^ (2.10) 

k a ±k b± 

and since it depends on k a ±kb_i, it is not a constant within the integral. However, the dominant, 
i.e. the leading logarithmic, contribution to Eq. ( |2.6| ) comes from the largest value of Ay, which 
is attained at the transverse momentum threshold, thus in Ref. |Q] Ay is fixed at its maximum 
by reconstructing the x's at the kinematic threshold for jet production and setting them in a 
one-to-one correspondence with the jet rapidities 

*r = ^=e* xr = ^ e ~ Vb - (2 - n) 
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Then the factorization formula fl2.2| ) is determined at fixed x^ N ,x^ N . Having fixed the rapidity 
interval ( |2.10| ) to 



Ay = ]n Xa \ b 



(2.12) 



the integration over k a ± and k b x can be straightforwardly performed^, and the gluon-gluon cross 
section ( |2.5[ ) becomes 



Vgg{k a ±>E ± ,k b± >E A 



^g^(A 9 ,l) 



with 



1 

2^ 



e u>{y)z 

dv ~o 777 

-oo + 1/4 



cos(2i/lnt) . 



For a s Ay <C 1, we can expand Eq. ( 231 ) and obtain 

&gg (K± > E ± , h± > Ex ) 



1X01011 [l + 0((a s Ay) 2 )} 



2Ei 



(2.13) 



(2.14) 



(2.15) 



On the other hand, for Ay > 1 we can perform a saddle-point evaluation of Eq. ( |2.14| ), and 
using the small-f expansion (|2.8f ), we obtain the asymptotic behaviour of the gluon-gluon cross 
section 10, 01 



^Hk a x>Ex,k b ±>Ex)- nCAas 6 



2Ej_ irBAy/4 



(2.16) 



At very large rapidities the resummed gluon-gluon cross section grows exponentially with Ay, 
in contrast to the LO (0(a 2 )) cross section ( |2.15| ), which is constant at large Ay. From the 
asymptotic formula (|2.16j ) the effective BFKL intercept a BFKL = A + 1 can be derived. In the 
experiment of Ref. [18] the BFKL intercept is measured by considering the ratio of hadronic 



cross sections, Eq. (2.2), obtained at different centre-of-mass energies and at fixed xip and scale. 
This allows the dependence on the parton densities to cancel, and the ratio of hadronic cross 
sections is therefore approximately equal to that of partonic cross sections evaluated at the 
relevant Ay values. 

In Eqs. ( [2.13| )-( ^.16 ) we have summarised the standard Mueller-Navelet analysis in which 



it is assumed that the x's are reconstructed through Eq. ( 2.11 ), and that the jet transverse 
momenta are unbounded from above. However, this is not the case for the DO analysis, since 

a) DO collect data with an upper bound on Q 2 = k a xk b x, which is of the same order of 
magnitude as the square of the lower cut on the jet transverse momenta, and thus cannot 
be ignored in the integration over the transverse momenta. 



b) DO reconstruct the x's through Eq. ( |1.2| ), which is well approximated by Eq. (2.3), but not 
by Eq. (EIH) ; 

We examine these two issues, and the modifications they entail on Eqs. ( |2.13| )-( ^.16| ), in turn. 
^In order to do the integrals analytically, it is necessary to fix the factorization scale fiF in Eq. ( |2.2[ ), e.g. 

IMF = E±. 
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2.2 Dijet production with an upper bound on Q 2 



In the Mueller-Navelet analysis the integration over the transverse momenta is taken up to 
infinity on the grounds that a finite and large upper bound on the transverse momenta would 
entail a contribution which is power suppressed in the ratio of the jet threshold to the upper 



bound. However, DO collect data with an upper bound on Q = k a ±kb±, namely Q 2 



1000 



GeV 2 , while Q 2 mia = Ej_ = 400 GeV 2 . When Q 2 ahl ~ Q 2 ^^, tne u PP er bound cannot be ignored 
in the integration over the transverse momenta. 

In order to assess what the modification on Eqs. ( 2.13] )~ (|2,16] ) is, we integrate the gluon- 
gluon cross section ( |2.5| ) over the transverse momenta k a ± and kb± above a threshold E± with 
the upper cut Q 2 < (Jmax imposed. We obtain 



&gg(ka± > E±, k b± > E±, k a± k b± < Q 



2 ) 

max/ 



nC 2 A a 2 
~2E 2 ~ 



^(Ay,l) 



Q 



* U * 



2 

max 



Q 



2 

max 



with rapidity interval Ay defined in Eq. (2.12), T defined in Eq. ([2.140 , and 

e u>(y)z S i n (2i/lnt) 



G(z,t) = F(z,t) 



1 

2^ 



dv 



i/ 2 + l/4 2v 



(2.17) 



(2.18) 



The analytic form of Eq. ( 2.17 ) depends on the particular definition of the upper cutoff Q max that 
DO uses, and changes substantially the shape of the gluon-gluon cross section (see Fig. ||) and 
in particular its subasymptotic dependence on Ay. At a s Ay <C 1, we expand the exponentials 
m Eq. Q2.17D and obtain 



cr. 



gg\ n a- 



> E ± , k b ± > E ± , k a± k b± < Q ri 



TlC 2 A 0? s 

~2ET 



E 2 



Q 



2 

max 



+ O (a s Ay) 



(2.19) 



Thus for the DO cuts, Q max = 1000 GeV 2 , Eq. ( gig ) lowers the LO cross section (|2~T5| ) by 40%. 
At Ay > 1 a saddle-point evaluation of Eq. ( p. 17 ) yields 



r(Ay>l), 

99 



' ij (k a ± > E ± , k b± > E ± , k a± k b± < Q 



2 ) 

max/ 



s~<2 2 
7TL 'A a s AAy 

2E 2 



(2.20) 



,BAy/4 



( HQ 2 ma jE 2 ± )+BAy/2 \ _ 
. V ^/BAy^ J 



Q 



2 

max 



V VBAy- 



with the error function 



3>(x) 



dte~ l . 



(2.21) 



Using the asymptotic expansion of the error function at x 3> 1, for very large rapidities ^JBAy 3> 
1 and at fixed Qmax/^'i we obtain 



a 



99 



(k a ± > E ± , k b ± > E ± , k a± k b± < QLix) 



irC\a 2 s 



D AAy 



2E{ y/TrBAy/4 



El 



Q 



2 

max 



o 2 

l + l n !lsp 
E 2 



(2.22) 

which is simply the asymptotic cross section of Eq. ( 2.16 ) reduced by a constant factor. For 
the DO values of E 2 ^ and Qmax' this corresponds to a reduction by a factor of about 4.3 in the 
standard asymptotic formula ( 2.16; ). 
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Although it might appear from Eq. (2.22) that the only effect of the Q max cut is to change 
the normalization relative to Eq. ( [2.16 ), which would drop out of the ratio of cross sections, one 
has to keep in mind that both equations are derived (from Eqs. fl2.17| ) and ( |2.13| ) respectively) in 
the asymptotic limit Ay 3> 1. In Fig. |l| we plot both the integral formulae and their asymptotic 
solutions. We see that the differences a gg — <5"g^ y>1 ^ are roughly constant with respect to Ay, 
and thus the relative differences get smaller with increasing Ay. However, at the Ay values 
relevant to DO analysis, it appears that non-negligible subleading corrections to the asymptotic 
formulae should be taken into account when determining the effective BFKL intercept. As can 
be inferred from Fig. |], these effects are more important when a Q 2 cut is imposed, since in this 
case it takes longer for the exponential rise with Ay to set in. 

In conclusion, the effect of an upper bound on the product of the jet transverse momenta can 
have a significant effect on both the normalization and Ay dependence of the gluon-gluon cross 
section. For the DO values, the increase of the cross section from small to large Ay is weakened 
by a factor of approximately 2, as shown in Fig. Care must therefore be taken in attributing 
any observed cross section increase exclusively to the e A ^ y 'BFKL' factor. 
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Cross section vs. rapidity difference 



~i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



+ BFKL 
x BFKL w. cuts 
— (2.17) 

(2.13) 

(2.22) 

(2.16) 




4 5 6 

Rapidity difference 



Figure 1: The dependence of the LL BFKL gluon-gluon subprocess cross section on the 
dijet rapidity separation Ay, without (Eq. ( 2. IS ), upper dashed line) and with (Eq. ( 2.l\ ), 
lower solid line) the Q^ax cu ^- The 'data points' are the same quantities calculated using 
the BFKL Monte Carlo discussed in Section with a s fixed and no additional kinematic 
cuts. Also shown are the asymptotic Ay 1 approximations, Eqs. (2. It) and ( 2.2\ ). 
The parameter values are a s = 0.164, E± = 20 GeV, Q t 2 nax = 1000 GeV 2 . 
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2.3 Dijet production at x's fixed as in the DO set-up 



In the analysis performed by DO, the x's are reconstructed through Eq. Ql.2| ). Since the jets are 
selected by requiring that Ay > 2, Eq. Q2.3] ) is a good approximation to Eq. Ql.2j). Conversely, 
the x's ( [2.11 ) used in the Mueller-Navelet analysis are by definition a good approximation to 
the DO x's ( |l.2[) only at the kinematic threshold for jet production. Therefore in this section we 
shall examine the modifications induced on Eqs. ( 2. 13| ) — ( 2.16| ) by defining the x's as in Eq. ([2.3]). 
First, we note that in this case the jet rapidities are not fixed, rather in a given x®) bin 
all the transverse momenta and rapidities contribute which fulfil Eq. (2.3). Thus the rapidity 
interval between the jets cannot be used as an independent, fixed observable. For convenience, 
we rewrite the rapidity interval ( |2.10| ) as 

E 2 



Ay = Y + In 



k a ±h± 



withQ 



Y 



ln^^ 
E 2 , 



(2.23) 



(2.24) 



The requirement that the rapidity interval be positive, Ay > 0, imposes an effective upper 
bound on Q 2 , 



Q 



2 

max 



2 n Y 



Eje 



Integrating then the gluon-gluon cross section (|2 
and fixed coupling a s , we obtain 



(2.25) 

over k a ± and k b ± above E±, at fixed x° a , x® 



v gg (ka± > E±, k b ± > E±, k a± k b ±< Q 
with 



2 ) 

max I 



TtC 2 A a 2 s 
~2El~ 



HY, 1) 



E 2 



Q 



2 

max 



G[Y, 



E] 



Q 



2 

max 



f(z,t) 



1 

2^ 



dv ■ 



V 2 + 



[l+u(u)Y 



t^ v) cos(2i/lnt), 



(2.26) 



(2.27) 



G(z,t) 



dv ■ 



V 2 + 



sin (2z/lnt) 



2v 



(2.28) 



and w{v) as in Eq. (2/7). Note that as Y — ► in Eq. ( 2.25 ), the upper bound on Q 2 goes to the 
kinematic threshold, Q max — » E"j_, and accordingly the cross section ( p.26[ ) vanishes. Note also 
that the tilde functions T ( p. 27 ) and Q (2.28) reduce to the functions J- ( p,14| ) and Q ( p. 18 ) for 
a s <C 1 but a s y ~ 1. This is understandable because the limit a s <C 1 is equivalent to neglecting 
subleading corrections, since (a s ) 2 Ay <C a s Ay, and therefore to neglecting differences in the 
definition of the rapidity interval, Ay ~ Y. For F>lwe perform a saddle-point evaluation of 
Eq. ( 2.26| ), and obtain the asymptotic behaviour 



99 



<i',:„ ^'(k a ±>E±,k b± >E_ L ,k a± k b± < Q 



2 ) 

max / 



irC 2 A a 2 



,AY 



2E{ ^/nBY/4 (1 + Af 



1 



E] 



Q 



2 

max 



1+A 



o 2 

1 + (1 + A) In ^5fS 



(2.29) 



*The constant Y resembles the rapidity interval ( |2 . 1 2| ) used in the Mueller-Navelet analysis, however it is not 
the same since Eq. (2.23) entails that Ay < Y , while Eqs. (2.1C) and (2.12) are just two different ways of defining 
the same rapidity interval. 
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We can also use the above analysis to include the DO experimental cuts of Q 2 < 1000 GeV 2 
and Ay > 2. In this case the analysis holds unchanged except that the upper bound on Q 2 is 
given by 

Q max = min(1000 GeV 2 , E\ e {Y ~ 2) ) , (2.30) 

where we have used the fact that Ay > 2 imposes the second effective upper bound on Q 2 . The 
shape of the cross section as a function of Y depends crucially on whether the upper bound on 
Q 2 is given by Eq. ( |2.25| ) or ( 2.3C| ) (see Fig. ^). This is more clearly apparent in the asymptotic 
region, F> 1, since for the upper bound ( |2.25| ) we can safely take (5 max — * 00 > with only the 
first term in the square brackets of Eq. (|2.29[) contributing; conversely, when the upper bound 
is given by Eq. (|273C| ), the sharp cutoff Q max = 1000 GeV 2 is much more restrictive than the 
bound ( 2.25| ) and depletes the cross section, which is given by the whole Eq. (|2729| ) . 

Using Fig. ||, we can get some idea of the expected effect of the Ay > 2 cut on the cross 
section ratio measured by DO. From Eq. ( |2.30| ) we see that this cut is inconsequential when 



Y > 2 + ln(1000 GeV 2 /E\) ~ 2.92 . 



(2.31) 



where we have used E± = 20 GeV. Conversely, this cut removes the entire cross section for 
Y < 2. For y/S = 1800 GeV we find Y > 2.92 for all bins, so the cut has no effect. However, for 
= 630 GeV we find Y = 2.37 in bin 1, Y = 2.63 in bin 2, and Y = 2.88 in bin 3, where we 
have used the minimum x\ and X2 in each bin to evaluate Y. Thus, bins 1 and 2 (and to some 
extent bin 3) get depleted at 630 GeV, simply due to the Ay > 2 cut. In section 4 we will see 
that this leads to a large cross section ratio in these bins, independent of the BFKL dynamics. 



Finally, we note that the asymptotic cross section Eq. ( 2.29| ) has the same shape in Y as 
Eq. (2.16) in Ay but different normalization: at a s (Q 2 = 400GeV 2 ) = 0.164, the normalization 
of Eq. ( jT29D with upper bound ( g25p is a factor 2.1 smaller than the one of the standard 
asymptotic formula ( p. 16 ), which becomes a factor 5.4 smaller than the one of Eq. (2.16) if the 
upper bound ( 2.30| ) is used. However, as Fig. ^ shows, for the values of rapidity interval relevant 
to the DO analysis we are far from the asymptotic region, and thus all the caveats made at the 



end of Section 2.2 on the extraction of the BFKL intercept from the DO data apply in this case 
as well. 



3 Dijet production and energy-momentum conservation 

Besides the effects described in the previous section, there are other facts that need to be 
considered when aiming at comparing experimental data to BFKL predictions, some of general 
nature and some specific to the experimental set-up of Ref. fTfj| . We will comment on the latter 
in the next section. As far as the former are concerned, we remind the reader of the following 
facts. (£) The LL BFKL resummation is performed at fixed coupling constant, and thus any 
variation in the scale at which a s is evaluated appears in the next-to-leading-logarithmic (NLL) 
terms, (ii) Because of the strong rapidity ordering, any two-parton invariant mass is large. Thus 
there are no collinear divergences in the LL BFKL resummation; jets are determined only at 
LO and accordingly have a trivial structure, (iii) Energy and longitudinal momentum are not 
conserved, and since the momentum fractions x of the incoming partons are reconstructed from 
the kinematic variables of the outgoing partons, the BFKL theory can severely underestimate 
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Dijet cross section vs. rapidity difference 




1 2 3 4 5 6 

A y or Y 



Figure 2: The dependence of the LL BFKL gluon-gluon cross section on Ay in the 
standard Mueller- Navelet calculation (Eq. ( 2.19Q ) (upper solid line) and on Y for the DO 
setup (Eq. \2.2<\ )). Four curves are shown for the definition of x 's applied in the DO 
analysis: Dashed line for the requirement Ay > 0, dotted line for Ay > 2, dash-dotted for 
Qmax of Eq. (2.3L ) and finally the lower, fat dash-dotted line for the asymptotic behaviour 
(Eq. (2.2$)) using Q„ lax of Eq. (2.3L ). The histograms are filled using the MC. 



the exact value of the x's, and thus grossly overestimate the parton luminosities. In fact, if n + 2 
partons are produced, energy- momentum conservation gives 



e Va / v-^ 



x b = ^ ( k bX + k aX e~ Ay + J2k iX e-y*+y» ] . (3.1) 



i=l 



The momentum fractions ( |2,3D in the high-energy limit are recovered by imposing the strong 
rapidity ordering ( |2.1[ ). However, the requirement x a ,Xb < 1 effectively imposes an upper limit 
on the transverse momentum (ki±) integrals. This effect is completely analogous to that we 



have considered in sections 2.2 and 2.3 



In an attempt to go beyond the analytic leading-logarithm BFKL results, a Monte Carlo 
approach has been adopted [25, 10, p6|| . If we transform the relevant Green's function to moment 
space via 

f(q a± ,q bx ,Ay) = J ^ Uq a± , q b± ) , (3.2) 
(by averaging / over azimuthal angles one obtains the quantity / introduced in Eq. (|2.5D) we 
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can write the BFKL equation as 

w fw(q a ±,qb±) = 

where the kernel K is given by 

K{q a ±,qb±, k±) = f w {q a L + k±,q b ± 



1 xr- - \ , Qs 
^o{q a± -96±) H 

2 7T 



cr/c i 



K(q a ±,qb±,k±) 



(3.3) 



(3.4) 



*1 + (<7a± + &±) 2 

The first term in the kernel accounts for the emission of a real gluon of transverse momentum 
k± and the second term accounts for the virtual radiative corrections. By solving the BFKL 
equation ( |3.3| ) by iteration, which amounts to 'unfolding' the summation over the intermediate 
radiated gluons and making their contributions explicit, it is possible to include the effects of 
both the running coupling and the overall kinematic constraints. It is also straightforward to 
implement the resulting iterated solution in an event generator. 

The first step in this procedure is to separate the k± integral in (|3.3j ) into 'resolved' and 
'unresolved' contributions, according to whether they lie above or below a small transverse 
energy scale fi. The scale fj, is assumed to be small compared to the other relevant scales in 
the problem (the minimum transverse momentum E±_ for example) . The virtual and unresolved 
contributions are then combined into a single, finite integral. The BFKL equation becomes 

d 2 k ± 



u f<j(qa±,ql±) 



1 xt-> - ^ , ° s 
7T 0{q a ± -qb_L) H 

Z 7T 



k± >fi 2 



k 2 . 



fu>{q a ± + k±,ql±) 



a s 
+ — 

7T 



d 2 k ± 



fcu(qa± + k±,q b ±) e(fj, 2 - k\) 



ql_Lfu{qa±,qb_L) 

k\ + (q a ± + k ± y 



. (3.5) 



The combined unresolved/virtual integral can be simplified by noting that since k\ <C q\±-,q\\_ 
by construction, the k± term in the argument of f u can be neglected, giving 



(w - u ) fu(q a ±,Qb±) 



1 



S(q a ± -ql±) + 



7T 



where 



u = — 

7T 



d 2 k A 



9(^ 2 - k 2 ± ) 



q 2 a ± 



d 2 k ± 



fUq a ±_ + k±,q bl _) 



a s In 



(3.6) 



(3.7) 



kl + (q a ± + k ± ) 2 _ 

The virtual and unresolved contributions are now contained in ujq and we are left with an 
integral over resolved real gluons. We can now solve ( |3.6| ) iteratively, and performing the inverse 
transform we have 



f(q a ±,qb±,&y) = 53/ (n) (go±,9fc±,Ay) ■ 



(3.8) 



n=0 



where 



f (0) (qa±,qb±,&y) 

f {n - 1] {q a L,qbL^y) 



k 2 ~ 


a s Ay 


{&.. 




k 2 ~ 


a s Ay 






Trk 2 , 


0(kh 



aX 



qb± 



n r -, I i n 

Y\ / d 2 k i± dyi Fi \ - 5(q a ± - q b ± - ^ fcjjj , 



.i=l 
,2 



fj, )0(yi-i - yi) 



(qa± + ET=A) 2 

(q a ± + E}=i kj±) 2 



1=1 

OisVi 



(3.9) 
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Thus the solution to the BFKL equation is recast in terms of phase space integrals for resolved 
gluon emissions, with form factors representing the net effect of unresolved and virtual emissions. 
In this way, each depends on the resolution parameter fi, whereas the full sum / does not. 
Changing simply shifts parts of the cross section from one to another .| 

Unlike in the case of DGLAP evolution, there is no strong ordering of the transverse momenta 
ki±. Strictly speaking, the derivation given above only applies for fixed coupling because we 
have left a s outside the integrals. The modifications necessary to account for a running coupling 
a s {kf 1 ) are straightforward [fTofl - 

The expression for / in Eqs. ([[!]) and (EO]) is amenable to numerical integration which can 



be implemented in the form of a Monte Carlo event generator |25| , |10| , 26 1. In this way we can 
trivially impose further cuts and combine the subprocess cross section with parton densities. 
Having made explicit the BFKL gluon emission phase space, overall energy and momentum 
conservation is imposed by using the momentum fraction definitions in Eq. (3.1). We can also 
run the Monte Carlo in 'LL mode', i.e. with fixed a s and no additional cuts or parton densities, 
to check that it does indeed reproduce the analytic results of the previous section, see Fig. |l[ 



4 Equal transverse momentum cuts: a dangerous choice 



We now take a closer look at the set-up specific to the DO analysis of Ref . [ 18 1 . As a preliminary 
observation, we might say that the values Ay probed are quite far from the asymptotic region 
where Eq. ( 2.16| ) is expected to hold, particularly at \/S = 630 GeV, where Ay is of the order 



of 2 to 3; unfortunately, here the only solution is to wait for the LHC to come into operation. A 
more serious, but solvable, problem is the following: dijet rates are quite sensitive to the emission 
of soft and collinear gluons, in the case in which they are defined by imposing equal cuts on 
the transverse energies of the two tagged jets. In this sense, a dijet total rate is completely 
analogous to the azimuthal correlation mentioned above. A detailed discussion on this point is 
given in Ref. , and will not be repeated here. In the current study, we will limit ourselves to 



illustrating the discussion of Ref. [27] by means of examples relevant to dijet production at the 
Tevatron. We will do this in two steps. First, in Subsection |4.1| , we will study this issue using 
a fixed-order perturbative computation, showing that dijet cross sections defined with unequal 
transverse momentum cuts do not have the same problems as those defined with equal cuts. 



Then, in Subsection 4.2, we will repeat the analysis of Section in the more general case of 



unequal transverse momentum cuts. 



^One other technical point deserves comment. The decomposition of / into the individual / (n) breaks the 
symmetry between q a ± and qt± which is only restored in the sum over n. Care must be taken when imposing 
asymmetric cuts on the a and 6 jets (see following section). Numerically, it is found that the fi dependence is 
slightly weaker when the higher E± cut is imposed on jet a. This asymmetry between the two jets is an artifact 
of the recursive solution. 



12 




Figure 3: Dijet rates, as defined in Eq. ( 4.1), for various cuts C. The cases of \/~S = 
630 GeV (left) and of y/S = 1800 GeV (right) are both considered. Dotted curves and 
circles have been rescaled by factor of 10 (left) and 50 (right). See the text for details. 



4.1 Dijet cross sections at fixed perturbative order 

In this subsection^], we consider jet production at fixed perturbative order in QCD. In particular, 
we use the partonic event generator of Ref. |^8| , which is accurate to NLO for any one- or two-jet 
observables. Similarly to what has been done previously for the gluon-gluon cross section, we 
define a total dijet cross section as follows: 

a(V,C) = a{k a± >E ± ,k b± >E ± +V,C), (4.1) 

where C generically indicates a set of cuts to be added to transverse-momentum cuts. As already 
mentioned, DO have 

C : \yi\ < 3, Ay > 2, Q 2 < 1000 GeV 2 (4.2) 

(i = a,b), together with some additional cuts on x\ and X2', furthermore, E± = 20 GeV. 

The rates defined in Eq. ( |4.ip are shown in Fig. ^, the left (right) panel presenting the case 
of yS = 630 GeV (y/S = 1800 GeV). Each plot consists of three sets of results, corresponding 
to different choices of C; for each of these choices, both the NLO results (displayed by the solid, 
dashed, and dotted curves) and the LO results (displayed by the boxes, diamonds and circles) 
are given. The solid curves and the boxes are obtained by imposing only the pseudorapidity cuts 
\yi\ < 3. The dashed curves and the diamonds correspond to the previous cuts on yi plus the 
cut Ay > 2. Finally, the dotted curves and the circles are relevant to the cuts given in Eq. (fh^), 
plus those that define bin number 1 (see Table ||). Notice that the results relevant to bin 1 have 
been multiplied by a factor of 10 and 50 at V~S = 630 GeV and V~S = 1800 GeV respectively, 
so that they can be shown together with the other results on the same plot. 

From the cross section definition in Eq. ( |4.1| ), it is clear that the smaller T>, the larger the 
phase space available; thus, one naively expects that the smaller T>, the larger the cross section. 

^In the fixed order pQCD analysis, at variance with Section [| we set fi F = {k a ± + ki,±)/2, a choice more suited 
to a Monte Carlo approach. 
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NLO 


Born 




V = 


V = 2 


V = 4 


V = 


V = 2 


V = 4 


bin 1 


-0.12(6) 


1.55(4) 


2.62(7) 


1.681(3) 


2.340(7) 


4.16(3) 


bin 2 


-0.16(4) 


1.14(3) 


1.46(4) 


1.265(3) 


1.417(4) 


1.739(8) 


bin 3 


-0.16(5) 


0.92(3) 


1.13(4) 


1.074(3) 


1.098(5) 


1.138(6) 


bin 4 


-0.19(6) 


0.92(4) 


1.15(5) 


1.036(4) 


1.045(6) 


1.068(8) 


bin 5 


-0.35(5) 


0.82(3) 


1.01(4) 


1.026(4) 


1.027(6) 


1.020(7) 


bin 6 


-0.45(9) 


0.82(7) 


1.08(9) 


1.015(9) 


1.01(1) 


1.00(1) 



Table 2: Fixed-order predictions for the ratio defined in Eq. (1.1). Numbers in paren- 
theses give the statistical error, which affect the last digit of the results shown. The values 
of T> are given in Ge V. 



This is indeed what happens at the LO level, regardless of the cuts C. On the other hand, 
the NLO cross section increases when T> decreases only if T> is not too close to zero; when T> 
approaches zero, o~{T>) has a local maximum, and then turns over, eventually dropping below the 
LO result. As discussed in Ref. at V = the NLO result is finite (i.e., does not diverge), but 
the slope da/dD is infinite. Fig. || thus clearly shows that at T> = (which corresponds to the 
definition adopted in the experimental analysis) the cross section is affected by large logarithms, 
that can spoil the analysis performed in terms of BFKL dynamics, exactly as in the case of the 
azimuthal decorrelation. 

Fig. H| already suggests a possible solution to this problem: simply define a dijet rate by 
considering different transverse momentum cuts on the two jets (that is, T> > 0). From the 
plots, we can expect that the potentially dangerous logarithms affecting the region T> = are 
not large starting from D of the order of 3 or 4 GeV. The figure might also at first sight seem 
to imply that a similar problem arises in the large T> region in the case in which the (physically 
relevant) cuts on x\ and xi are imposed (dotted curves and circles). However, it is easy to 
understand that in such a case the large difference between the NLO and LO results is simply 
due to phase space: in fact, at LO T> > effectively forces both jets to have k± > E± + T>; at 
NLO, this is no longer true. 

Let us therefore consider again the ratio of Eq. ( |1.1| ) , now rewritten to indicate explicitly the 
cuts adopted: 

R(V, C) = a{V, C;VS = 1800 GeV)/a(V, C;VS = 630 GeV) , (4.3) 

with C given in Eq. (f4.2j), and additional (binning) cuts on x\ and X2- Our predictions for R, 
both at NLO and LO, are presented in Table ^, where we show the results for all of the bins 
of Table |l[ The entries relevant to T> = display a pathological (negative) behaviour at NLO. 
However, even if they were positive, they could not be considered reliable, since any fixed-order 
QCD computation (beyond LO) is unable to give a sound prediction in this case. On the other 
hand, we see that for larger values of T> the situation improves, in the sense that it reproduces 
our naive expectation: the ratio should converge towards one, for increasing Ay (i.e., larger bin 
numbers); while for T> = 2 GeV the NLO results are still sizeably different from the LO results, 
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= 1.8 TeV 


y/S = 0.63 TeV 




V = 


V = 2 


V = 4 


V = 


V = 2 


V = 4 


bin 1 


31.01(3) 


21.24(3) 


14.19(2) 


18.46(3) 


9.080(2) 


3.399(1) 


bin 2 


22.66(2) 


15.52(2) 


10.36(2) 


17.91(2) 


10.98(2) 


5.969(2) 


bin 3 


13.50(2) 


9.22(2) 


6.16(2) 


12.57(2) 


8.39(2) 


5.41(2) 


bin 4 


12.11(3) 


8.29(2) 


5.54(2) 


11.69(3) 


7.20(2) 


5.19(2) 


bin 5 


7.19(1) 


4.90(1) 


3.26(1) 


7.10(1) 


4.78(2) 


3.19(1) 


bin 6 


4.25(2) 


2.89(2) 


1.92(1) 


4.19(2) 


2.85(2) 


1.92(1) 



Table 3: Cross sections in nanobarns as given in Eq. ({.1), at the LO and for two 
different centre- of -mass energies. Statistical errors are given in parentheses. 



in the case of T> = 4 GeV the NLO and LO results are statistically compatible (within one 
standard deviation) for bins 3-6, and they are both approaching one. 

Inspection of Fig. |3| and Table [2] tells us that, in order to avoid the presence of large logarithms 
of non-BFKL nature in the cross section, a value of T> = 4 GeV is probably a better choice than 
T> = 2 GeV. Of course, the larger T>, the smaller the cross section, and therefore the fewer the 
events. In order to give an estimate of the loss of events that one faces when going from D = 



to larger values, we present in Table ^ our LO predictions for the rate defined in Eq. (4.1), with 
the cuts of Eq. d4.2|) and our six bins. Of course, it is well known that NLO corrections are 
mandatory in jet physics to get good agreement with data. However, here we just want to have 
a rough idea of the number of events lost when increasing one of the transverse energy cuts; this 
number is sensibly predicted by the ratio a{T>)/a(T> = 0), even if a is only computed at LO. 
From the table, we see that at T> = 4 GeV the number of events decreases, compared to the 
case T> = 0, by a factor slightly larger than two; this factor gets much larger only for the first 
two bins at yS = 630 GeV, which are however less relevant from the point of view of BFKL 
dynamics. 

4.2 Dijet production in the BFKL theory with an asymmetric cut 

We now turn again to the BFKL equation, and study the dependence on the offset T> introduced 
in the previous subsection, in essentially the same way as we did in Section ||. We start by 
integrating the gluon-gluon cross section ( |2.5| ) over k a ± > E± and k b ± > E± + V with the upper 
cut Q 2 = k a ±k b ±< 1000 GeV 2 , and with the x's defined as in the Mueller-Navelet analysis, 
Eq. 



7rC 2 G? 

cr gg (k a± > Ej_, k b± >E L +V, k a± k b± < QLx) = 2E±{e"+ V) (44) 



T Ay, 



E± \ E ± (E ± +V) 



E ± + Vj 2QI 



E l\ n ( \ (E± + V) 2 



giAy,—^\ + giA y 



Ql ' ^ I — »' Q2 

Vinax/ \ MJn 



with T and Q defined in Eq. ( |2.14 ) and Eq. ( [2.18 ). Repeating the calculation as in Section 2.3 



with the x's defined as in Eq. (p.3j), yields a cross section of the same form as Eq. ( |4.4[) up to 
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replacing the rapidity interval ( |2.12| ) with the constant ( 2.24| ), the upper bound Qmax above 
with Eq. ( ^30| ) and the function T (G) with T (§), Eq. ( gg7| ) (Eq. (gjg)). At £> = we recover 
Eqs. ( p. 17 ) and fl2.26j ) respectively. However, near T> = Eq. (^4|) and its analogous one with 
the tilda functions display the same qualitative behaviour when expanded to NLO as the exact 
NLO cross section [p7|1 . In order to see this, we take Eq. ( fP[ ) in the limit Qmax — * °°> such that 
only the first term on the right hand side of Eq. ( |4.4[) survives. We analyse its NLO term by 
expanding its exponential to 0(a s ), 



Vgg(k a ± > E±, k b± >E ± +V) 



(4.5) 



C\al 



00 1 - 2a s Ay [Re^ (1/2 + iv) - V>(1)] 



2iu 



4E ± (E ± +V)J_ ao '~ ^ + 1/4 V^+W +0 «^) S 

The denominator has poles at v = ±i/2. For the LO term, the integration over is straightfor- 
ward. For the NLO term, we use the integral representation of the digamma function, 



^)-V(i) 



„z-l 



dx 



1 



and after performing the integrals over v and x, we find 

°gg(k a ± > Ej_, k b± > E ± + V) 



El 



1 - a s Ay 



-2E L V - V 1 

JeZ+W 



In- 



-2E ± V - V 2 



+ 2 In 



E ± +V 
Ex 



(4.6) 

(4.7) 
6{-V) 



+ 



. _ . (2E ± V + V 2 2E ± V + V 2 E ± 

1 — «,Au ^ m — -7T + 2 in 

y V E 2 (E ± + V) 2 E ± + V 



9(V) 



(E ± + V) 2 
+ O ((a s Ay) 2 ) . 

At LO, transverse-momentum conservation forces the cross section to behave like 1/-E^? for T> < 
and like 1/(E± +T>) 2 for T> > 0, even though the cuts over the transverse momenta are asym- 
metric, see Fig. |. For V = 0, Eq. (fO]) reduces to Eq. ( gig ). For |P| < Eq. ( p7| ) 
becomes 



1 



2\V\ 2\V\ , ' 
a s Ay-^ln-^ + 0(V) 

Il/J_ _C/_|_ 



+0 {(a s Ayf) . (4. 



The slope of Eq. (4.7) with respect to T> is negative (positive) for T> < 0(> 0), and infinite at 
T> = 0, in agreement with Ref. |p7| ]. In addition, by using Eq. (|4.8| ) to evaluate the ratio (|l.l| ) and 
remembering that asymptotically Ay A — ► Ay s + 2 ln(1800/630), we find that the NLO BFKL 
ratio also goes to 1 as Ay grows, in agreement with Table |2[ 

In the BFKL Monte Carlo approach, the implementation of asymmetric cuts on the jets is 
straightforward. For fixed a s and no additional cuts or parton densities, the analytic result of 



Eq. (4^4) is reproduced, see Fig. [|. 

Finally, we can use the BFKL Monte Carlo to calculate the 'DO' cross section ratios defined 
m Eq. (O), in the various bins. Table gives the predictions using the Monte Carlo run in two 
modes: 

Naive: fixed a s , no kinematic constraints, parton densities evaluated at Bjorken x's given in 
Eq. Q. 
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Naive BFKL 


BFKL MC 




V = 


V = 2 


V = 4 


V = 


V = 2 


V = 4 


bin 1 


2.16(4) 


2.56(5) 


3.61(8) 


1.615(3) 


2.013(5) 


2.922(8) 


bin 2 


1.47(4) 


1.53(4) 


1.73(4) 


1.048(2) 


1.114(2) 


1.289(3) 


bin 3 


1.22(4) 


1.16(4) 


1.10(3) 


0.866(2) 


0.851(2) 


0.872(2) 


bin 4 


1.18(4) 


1.12(4) 


1.14(4) 


0.806(4) 


0.783(2) 


0.787(2) 


bin 5 


1.26(6) 


1.14(5) 


1.11(5) 


0.847(2) 


0.824(2) 


0.820(3) 


bin 6 


1.31(8) 


1.28(7) 


1.20(6) 


0.863(3) 


0.841(3) 


0.838(3) 



Table 4: BFKL predictions for the ratio defined in Eq. (1.1). Numbers in parentheses 
give the statistical error, which affect the last digit of the results shown. The values of T> 
are given in GeV. 



Full: running a s , energy-momentum conservation applied, parton densities evaluated at x a , xj, 
values given in Eq. (|3.1|). 

Evidently neither the naive BFKL nor the BFKL MC calculation shows the 'pathological' be- 
haviour of the exact NLO calculation at T> = (this is already apparent from Fig. |j). Instead the 
numbers are quite stable against variations in T>. For all T> the naive BFKL calculation shows 
an initial decrease in the cross section ratios, before reaching a minimum around bin 4 where the 
expected rise due to BFKL dynamics sets in. The initial decrease is simply the subasymptotic 
effect of the Ay > 2 cut on the cross section at yfs = 630 GeV, as discussed in Section |2.3| , 
and consistent with the qualitative behaviour of the Born cross section ratios in Table ^[ On 
the other hand, if the effects from the parton densities did factorize out completely, one would 
expect an asymptotic (in bin number) ratio of R = (Sa/Sb) X with A = ""^^M ln 2 ~ 0.45 [26|. 



This gives R ~ 2.6 for the DO values. However, we have already argued that such a rise is 
not expected in the DO analysis, mainly because of the rather stringent Q 2 cut (see Fig. |] and 
Fig. ^). Furthermore, the expected rise is also slightly decreased by the cut in x a ,Xb < 1, and 
hence on k a ± and kbx, introduced by the parton densities. The full BFKL MC calculation ratios 
also show an initial decrease to a minimum around bin 4. However now the ratio is below 1 
already from bin 3 onwards. Such an effect was already reported and explained in Ref. pqj . It is 
a kinematic effect due to an effective upper limit on the transverse momentum allowed for each 
emitted gluon. As the rapidity separation between the dijets is increased towards its maximum 
allowed value, the BFKL gluon phase space is squeezed from above and the 'naive' cross section 
is heavily suppressed. The higher the collision energy the more dramatic the effect, and hence 
the ratio R falls below 1. 



5 Conclusions 

In this paper, we have reconsidered the suggestion by Mueller and Navelet of studying dijet 
cross sections at large rapidity intervals and for different hadronic centre-of-mass energies, in 
order to find evidence of BFKL physics. We were motivated by a recent paper by the DO 
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Figure 4: The dependence of the gluon-gluon subprocess cross section on the offset T>, 
for fixed separation Ay = 3. The resummed prediction Eq. (4-4) * s shown as a solid line, 
with the results of the corresponding BFKL Monte Carlo calculation superimposed. The 
dash- dotted line is the LO contribution, and the dashed line is the 0(a s ) contribution of 
Eq. (4/1). The parameter values are a s 



0.164, E x = 20 GeV, Q 



2 

max 



OO. 



Collaboration [18|, where dijet data have been used to measure the effective BFKL intercept by 
comparison with the standard analytic asymptotic formulae given by Mueller and Navelet. In 
fact, we have argued that the definition of the momentum fractions used by DO and some of 
the acceptance cuts imposed by DO spoil the correctness of this procedure, and require a more 
careful theoretical investigation. 

In particular, we are concerned by a difference between DO and the standard Mueller-Navelet 
analysis in the reconstruction of the momentum fraction x of the incoming partons, by the 
presence of an upper bound on the momentum transfer Q 2 , and by the requirement that the 
two tagged jets have the same minimum transverse energy. 

The Q 2 cut allows, at the experimental level, and together with the binning cuts on xi,2, a 
reduction of the systematic errors, since in the ratio of the cross sections measured at different 
centre-of-mass energies the dependence on the parton densities cancels to a significant extent. 
We have shown that, at the level of partonic cross sections, the upper bound on Q 2 and the x's 
used in the DO analysis reduce the Mueller-Navelet cross section by a factor of more than 5. On 
the other hand, the dependence on such a cut, as well as the dependence on the precise definition 
of the x's, cancel out when considering the ratio of cross sections obtained at different energies. 
However, this is only true when the asymptotic forms of the cross sections are considered. 
Unfortunately, at the energies and rapidity intervals probed at the Tevatron, it appears that 
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the asymptotic expansions do not reproduce accurately enough the exact analytic results; in 
particular, the quality of the approximations are worse in the case in which an upper cut on Q 2 
is imposed. We are therefore led to conclude that, regardless of the use of cross sections or of 
rates of cross sections to study BFKL physics, the effect of an upper bound on Q 2 cannot be 
ignored. 

As far as the cuts on the transverse momenta of the trigger jets are concerned, we have 
pointed out that in the case in which such cuts are chosen to be equal, the cross sections are 
plagued with large logarithms of perturbative, non-BFKL origin. In this sense, the total dijet 
rates are therefore on the same footing as the azimuthal decorrelations. We therefore believe 
that a much safer choice is to have different cuts on the transverse momenta of the two jets. 
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